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ABSTRACT 


A  brief  review  of  Task  II  of  the  research  program  on  "Three- 
Dimensional  Internal  Flows  in  Turbomachinery,"  at  the  University 
of  Cincinnati,  is  presented.  The  review  covers  the  period 
between  January  1979,  and  December,  1978. 

During  this  period,  five  related  areas  have  been  under 
investigation.  These  are  (1)  Flow  through  Joukowski  Ducts; 

(2)  Flow  through  Curved  Polar  Ducts;  (3)  Analysis  of  Streamwise 
Flow  Separation  via  a  Semi-Elliptic  Mathematical  Model; 

(4)  Higher  Accuracy  of  Numerical  Solutions  and  (5)  Increasing 
Efficiency  of  Numerical  Algorithms. 

The  personnel  involved  in  this  research  program,  the 
progress  during  the  year  as  well  as  further  work  planned  for 
the  remainder  of  the  contract  period  have  been  briefly 


outlined . 
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1.  REVIEW  OF  RESEARCH  PROJECT 

The  primary  objective  of  this  two-year  research  project  is 
to  develop  a  finite-difference  numerical  technique  for  the 
solution  of  the  parabolized  Navier-Stokes  equations  using  a 
system  of  surface-oriented  coordinates  suitable  for  ducts  of 
arbitrary  geometry.  This  objective  is  to  be  achieved  within 
the  framework  of  a  laminar-flow  analysis.  An  additional 
objective  to  be  fulfilled  by  this  research  consists  of  developing 
an  analysis  for  laminar  duct  flow  with  streamwise  separation 
using  mathematical  models  simpler  than  the  full  Navier-Stokes 
equations.  Finally,  a  finite-difference  numerical  technique  is 
to  be  developed  using  a  two-equation  turbulence  model  for 
analyzing  flows  in  straight  and  curved  ducts  of  regular  cross 
sections,  with  particular  emphasis  on  the  accuracy  of  the 
analysis  in  the  wall  region. 

In  order  to  achieve  the  above  objectives,  the  following 
three  problems  have  been  studied  during  the  first  year  of  the 
twe-vear  contract  period: 

i.  Laminar  flow  through  Joukowski  ducts, 
ii.  Laminar  flow  through  curved  polar  ducts, 
and  iii.  Analysis  of  laminar  streamwise  flow  separation  via 
a  semi-elliptic  mathematical  model. 

In  each  phase  of  the  research,  attention  is  given  to  the  accuracy 
of  the  numerical  solutions  as  well  as  the  efficiency  with  which 
they  can  be  determined.  For  this  purpose,  the  following  two 
additional  areas  have  been  studied: 
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iv.  Higher  accuracy  of  numerical  solutions, 


and  v.  Efficient  numerical  algorithms. 

Suitable  model  problems  have  been  selected,  wherever  required, 
in  this  research,  with  the  overall  goal  of  contributing  to 
further  the  understanding  of  internal  flows  in  turbomachinery. 

Each  of  the  above  research  areas  is  described  in  the  following 
subsections . 

1 . 1  Laminar  Flow  Through  Joukowski  Ducts 

The  flow  in  the  blade  passages  of  a  turbomachine  (compressor 
or  turbine)  is  three-dimensional,  unsteady,  viscous  and  turbulent. 
Further,  the  blade  passages  are  geometrically  complex.  A  thorough 
understanding  of  this  highly  complex  flow  is  needed  in  order  to 
achieve  an  optimum  blade  design  and  thereby  gain  improvement 
in  engine  performance.  Theoretical  techniques  that  can  fulfill 
these  objectives  can  serve  to  complement  and  partly  replace  the 
costly  experimental  triai-and-error  design  procedures.  Therefore, 
it  is  necessary  to  develop  theoretical  and  numerical  techniques 
for  analyzing  three-dimensional,  viscous  flow  in  non-simple 
geometries.  The  purpose  of  the  present  study  is  to  develop  an 
analysis  and  solution  technique  for  the  three-dimensional  viscous 
flow  in  ducts  of  arbitrarily  varying  cross  sections,  therefore, 
the  study  is  applicable  to  cascade-like  passages. 

The  geometry  of  the  cascade  channel  is  transformed  from  its 
complicated  shape  in  the  physical  domain  to  a  simple  rectangular 
computational  domain  where  the  flew  boundaries  align  conveniently 
with  the  transformed  coordinates.  The  coordinates  are  adjusted 


such  that  one  coordinate  direction  is  approximately  aligned 
with  the  predominant  direction  of  the  flow;  this  later  permits 
parabolization  of  the  transformed  flow  equations.  Because  of 
the  complex  nature  of  the  geometry,  the  transformation  between 
the  physical  coordinates  and  the  computational  coordinates  is 
established  through  a  numerical  procedure  (Ref.  1) .  Accordingly, 
the  transformed  coordinates  are  obtained  as  solutions  of  Poisson 
equations  in  the  physical  domain,  with  Dirichlet  boundary 
conditions.  Ever,  the  coordinate-transformation  eauations  are 
solved  in  the  transformed  domain  so  as  to  allow  the  use  of  a 
uniform  computational  grid.  For  this  purpose,  the  role  of  the 
dependent  and  the  independent  variables  is  interchanged  and  the 
inverted-transformation  eauations  are  then  solved  using  an 
alternating-direction  implicit  (ADI)  numerical  method.  The 
coordinates  so  generated  are  not  necessarily  orthogonal.  In 
the  present  study,  only  radially  invariant  geometries  have  been 
considered.  A  possible  system  of  transformed  coordinates  is 
shown  in  Fig.  1  for  a  typical  turbine  cascade. 

The  mathematical  model  of  the  flow  problem  is  derived 
from  the  Cartesian-  or  cylindrical-coordinate  continuity  and 
momentum  equations  transformed  to  the  general  non-orthogonal 
surface-oriented  coordinates  for  the  duct.  The  flow  problem  is 
assumed  to  be  parabolic  in  the  main-flow  direction  because  of 
minimal  streamwise  viscous  diffusion.  This  permits  the  flow 
equations  to  be  parabolized  with  respect  to  the  streamwise 
coordinate.  These  equations  then  lend  themselves  to  solution 
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by  a  marching  numerical  procedure,  thereby  avoiding  the  large 
computer-storage  and  time  requirements  for  solution  of  a  full 
three-dimensional  elliptic  system  of  equations. 

The  numerical  method  employed  for  obtaining  the  solutions 
consists  of  the  ADI  scheme  in  the  cross-plane,  with  a  Crank- 
Nicolson  type  marching  in  the  main-flow  direction.  Particular 
care  is  taken  to  preserve  conservation  of  the  mass-flow  rate 
across  the  duct  even  when  the  transformed  cross-sectional 
surface  is  not  completely  normal  to  the  duct  axis.  The  Neumann 
boundary-value  problems  encountered  in  the  analysis  are  solved 
numerically  using  the  SOR  (Successive  Over-Relaxation)  method 
by  points.  These  Neumann  problems  have  been  carefully  formu¬ 
lated  in  the  transformed  coordinates  so  as  to  obtain  a  rapidly 
converging  solution. 

Neither  experimental  nor  theoretical  results  are  available 
in  the  literature  for  3-D  laminar  viscous  flow  through  varying- 
area  ducts.  Hence,  in  order  to  assess  the  validity  of  the  present 
analysis  and  computational  procedures,  the  incompressible  viscous 
laminar  flow  through  a  two-dimensional  diverging  channel  is 
used  as  a  test  case.  An  exact  solution  of  the  Navier-Stokes 
eauations  for  this  problem  is  available  (Ref.  2)  in  a  polar 
coordinate  system.  The  results  obtained  by  the  present  apcroach 
show  very  good  agreement  with  the  exact  solutions. 

The  analysis  is  then  used  to  predict  the  three-dimensional 
viscous  flow  in  Joukowski  ducts  with  rectangular  as  well  as 
polar  cress  sections  (Fig.  2).  This  term  has  been  adopted  to 
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define  a  duct  of  rectangular  or  polar  cross-section  whose  aspect 
ratio  and  area  vary  with  streamwise  distance  due  to  the  presence 
of  a  Joukowski  airfoil  profile  on  two  of  its  opposite  walls.  The 
results  are  presented  in  terms  of  axial  variation  of  primary 
velocity,  mean-viscous-pressure  drop  and  cross-flow  velocity 
profiles.  The  results  obtained  are  also  analyzed  for  the  various 
problem  parameters  such  as  the  aspect  ratio,  Reynolds  number 
based  on  inlet  hydraulic  diameter  and  the  area  ratio.  Some 
sample  results  are  presented  in  Fias.  3-5.  In  all  the  cases 
studied,  the  constant-area  duct  beyond  the  airfoil  section  is 
extended  until  the  fully  developed  velocity  profile  is  obtained. 
The  case  of  a  polar  Joukowski  duct  is  a  very  good  approximation 
of  an  actual  turbine-stator  blade  passage.  For  a  typical  polar 
configuration,  the  development  of  axial  velocity  profiles,  the 
axial  variation  of  the  centerline  and  the  maximum  axial  velocities 
and  the  mean-viscous  pressure  drop  are  shown  in  Fig.  6. 

The  results  for  both  rectangular  and  polar  Joukowski  ducts 
are  found  to  be  in  conformity  with  the  general  physical  features 
of  the  flow  and  serve  as  a  strong  indication  of  the  adecuacy 
of  the  present  analysis  of  flow  in  geometrically  complex  con¬ 
figurations.  For  example,  in  the  vicinity  of  the  leading  edae 
of  the  airfoil,  the  presence  of  strong  secondary  flow  causes 
the  streamwise  velocity  to  attain  its  maximum  value  off  the 
centerline.  Consequently ,  between  the  ieadino-edge  section  and 
the  minimum-area  section,  the  axial-velocity  profile  shews  two 
peaks  m  it.  Also,  as  expected,  the  Reynolds  number  is  no 
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longer  simply  a  scaling  parameter  for  the  streamwise  distance 
even  for  the  axial-velocity  development.  Nevertheless,  the 
axial  velocity  as  well  as  the  mean  viscous-pressure  drop  are 
weak  functions  of  Re. 

A  technical  abstract  (Ref.  3)  summarizing  the  computational 
details  of  this  study  is  submitted  for  possible  presentation 
at  the  Second  International  Conference  on  Mathematical  Modelling 
in  July  1979,  in  St.  Louis,  Missouri;  a  written  version  of  the 
paper  will  then  appear  later  in  the  Conference  Proceedings.  A 
technical  abstract  (Ref.  4)  containing  the  main  body  of  the 
results  of  the  physical  flow  problem  is  submitted  for  presentation 
at  the  ASMt  1979  Winter  Annual  Meeting  to  be  held  in  New  York 
in  December  1979;  a  written  version  of  that  paper  is  to  appear 
in  the  Proceedings  of  an  ASME  Symposium.  Yet  another  pacer 
(Ref.  5)  based  on  this  work  is  presently  in  preparation  and  is 
to  be  submitted  for  possible  publication  in  a  refereed  technical 
journal . 

In  summary,  it  may  be  therefore  stated  that  the  study 
provides  a  viable  approach  for  the  accurate  prediction  of  the 
three-dimensional  viscous  flew  field,  including  the  recirculating 
secondary  flows,  in  ducts  with  arbitrarily  varying  cross-sections. 
Although  only  the  laminar  incompressible  case  has  been  cresentiv 
considered,  the  analysis  serves  as  a  significant  basic  step 
towards  the  analysis  of  three-dimensional  viscous  flew  fields  in 
turbomachinery  cascades.  For  applicability  to  a  real  cascade, 
several  extensions  and  generalizations  need  to  be  considered. 


These  include  additional  compl ications  in  the  geometry  such  as 
turning,  stagger  and  twis*  for  the  blades,  as  well  as  the 
implementation  of  the  periodicity  condition.  Additional  effects 
associated  with  the  physical  features  of  the  flow  are  those  due 
to  turbulence,  compressibility  and  upstream  influence  and  stream- 
wise  flow  separation.  All  these  are  important  effects  which 
can  now  be  considered,  one  by  one,  with  the  help  of  the  basic 
building  block  developed  in  the  present  work. 

1 . 2  Laminar  Flow  Through  Curved  Polar  Ducts 

A  large  and  important  class  of  engineering  problems  is 
represented  by  internal  viscous  flows  in  curved  ducts.  Accurate 
and  reliable  calculation  of  such  flows  is  of  great  practical 
interest  in  the  prediction  of  flow  through  inlets,  nozzles, 
diffusers  and  other  components  of  turbomachinery.  Curved  ducts 
of  polar  cross  section  (i.e.,  finite-annular  cross  section; 

Fig.  7)  are  of  particular  significance  in  many  of  these  appli¬ 
cations.  When  viscous  fluid  flows  through  the  curved  ducts, 
centrifugal  forces  induce  a  secondary  flow  in  the  crcss-sectional 
planes  and  lead  to  higher  pressure  losses  than  those  encountered 
in  ducts  without  longitudinal  curvature.  Thus,  curved-duct 
flows  are  truly  three-dimensional,  so  that  their  determination 
requires  the  solution  of  the  three-dimensional  elliptic  Navier- 
Stokes  equations,  while  analytical  solutions  for  these  flows 
are  impossible,  even  their  numerical  solutions  can  become  very 
expensive  in  terms  of  computer  costs.  Therefore,  other  methods 
are  sought  which  involve  some  simplification  of  the  Navier-Stokes 


equations,  with  minimum  compromise  on  the  class  of  problems 
that  can  be  treated.  Numerical  solutions  of  these  equations 
can  be  obtained  with  great  efficiency. 

In  the  present  study,  the  complex  three-dimensional  flow 
field  occurring  in  curved  polar  ducts  is  analyzed  and  determined 
numerically.  The  flow  problem  is  modelled  mathematically  in 
the  toroidal  coordinate  system  (r,9,o),  as  shown  in  Fig.  7, 
using  the  three-dimer. sional  Navier-Stokes  equations  parabolized 
in  the  longitudinal  direction  $.  Thus,  streamwise  diffusion 
terms  are  neglected  as  compared  to  the  cross-plane  diffusion 
terms.  This  also  leads  to  decoupling  of  the  axial  pressure 
gradient  from  the  cross-plane  pressure  gradients.  Parabolization 
requires  that  the  flow  Reynolds  number  is  not  small  and  that 
the  streamwise  velocity  is  nowhere  negative.  These  conditions 
are  well  satisfied  in  the  present  problem  provided  the  longi¬ 
tudinal  curvature  is  not  too  sharp.  This  formulation  permits 
a  marching  calculation  procedure  starting  from  a  given  set  of 
initial  conditions  to  the  solution  at  desired  downstream  locations. 

For  adequate  representation  of  the  high-gradient  regions 
near  the  duct  walls,  analytical  coordinate  transformations  are 
used  to  introduce  a  nonuniform  grid,  thus  circumventing  the 
need  for  excessive  computational  time  and  storace.  A  trans¬ 
formation,  which  had  been  successfully  used  earlier  for  the 
two-dimensional  cavity-flow  problem  (Ref.  6),  has  been  employed 
in  the  present  duct-flow  problem.  The  parameters  in  this  trans¬ 
formation  provide  control  over  the  first  step  near  the  wall. 

Thus,  the  first  step  can  be  varied  from  1/(2N)  to  1/(9N),  where  N 
is  the  total  number  of  intervals  employed. 
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The  solution  procedure  employed  parallels  that  used  in 
Refs.  7-9.  The  momentum  equations  are  solved  by  the  Douglas- 
Gunn  two-step  alternating-direction  implicit  method,  while  the 
Neumann  boundary-va lue  problems  encountered  in  the  analysis 
are  solved  using  the  successive  over-relaxation  method.  Since 
the  Poisson  operator  in  terms  of  the  present  transformed 
coordinates  has  varying  coefficients,  the  finite-difference 
representation  of  the  Neumann  problem  required  an  extreme  amount 
of  care  in  order  to  obtain  a  converged  solution.  Accurate 
satisfaction  of  an  integral  relation  between  the  source  term 
and  the  boundary  value  of  this  Neumann  problem  was  crucial  to 
obtaining  a  solution  of  the  problem. 

First,  the  reliability  of  the  present  analysis  and  the 
associated  computer  program  was  checked  by  reproducing  the 
results  for  some  cases  for  which  solutions  were  already  available. 
This  also  allowed  assessment  of  the  coordinate  transformation. 
Figure  3  3hows  typical  results  for  flow  in  a  curved  polar  duct  and 
a  curved  rectangular  duct  of  unity  aspect  ratio,  for  Dean  number 
K  *  55  and  100.  A  comparison  has  been  made  for  the  stream, wise 
variation  of  the  axial  centerline  velocity  and  the  mean-pressure 
drop  for  these  ducts.  Both  conf iguraticns  show  identical  behavior 
for  these  quantities.  Similar  conclusions  were  observed  earlier 
for  the  cor resoond mg  straight  polar  and  straight  rectanauiar 
ducts . 

Results  have  been  obtained  for  several  configurations  of 
the  curved  polar  duct,  with  various  values  of  the  problem  para¬ 
meters,  namely.-  Reynolds  number  Re,  Dean  number  K,  and  aspect 
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ratio  AR.  Figure  9  shows  the  streamwise  variation  of  the  skin- 


friction  coefficient  Cf  at  the  two  positions  denoted  by  A  and  B 

*w 

for  the  cross  section  shown.  For  each  case  presented,  C(  has 

its  maximum  value  near  the  initial  station  with  uniform  entrance 

flow  velocity  profile.  With  increasing  streamwise  distance, 

Cf  decreases  until  it  achieves  its  fully-developed  value  and 
rw 

remains  constant  further  downstream.  Also,  C*  decreases  with 

‘w 

increasing  Dean  number.  In  view  of  the  nature  of  centrifugal 

forces,  the  position  of  B  on  the  outer  surface  with  the  larger 

radius  of  longitudinal  curvature  and  the  position  of  A  on  the 

inner  surface  with  the  smaller  radius  of  transverse  curvature, 

Cf  at  B  is  always  larger  than  at  A. 

‘w  ‘w 

The  streamwise  development  of  the  axial  velocity  prefiles 

along  the  mid-surfaces  of  the  polar-duct  cross  section  is  shown 

in  Fig.  10.  It  should  be  observed  that  no  symmetry  exists  across 

either  of  these  mid-surfaces  in  the  curved  polar  duct.  The 

asymmetry  with  respect  to  6  is  due  to  the  longitudinal  curvature 

of  the  duct  and  is,  therefore,  significantly  more  pronounced 

than  the  asymmetry  with  respect  to  r  which  is  mainly  due  to 

transverse  curvature  effects.  As  expected,  the  degree  of  asymmetry 

increases  with  increase  in  the  Dean  number  K  as  seen  from  Fig.  ICa 

for  K  *  100  and  Fig.  10b  for  K  *  630. 

Typical  secondary-flow  velocity  profiles  are  shown  in 

Fig.  11  for  the  fully  developed  flow  in  two  curved  polar  duct 

conf igurations .  At  K  a  130,  Fig.  lia  shows  a  considerable  degree 

of  asymmetry  with  respect  to  this  is  further  enhanced  at 


K  *  630  (Fig.  lib).  Moreover,  at  K  *  630,  the  secondary  velocity 
profiles  indicate  additional  changes  in  secondary-flow  direction 
as  compared  to  the  case  for  K  *  100. 

The  results  indicate  that  the  study  provides  a  valid  analysi 
and  a  numerical  method  of  solution  for  flows  in  curved  ducts 
of  great  practical  significance.  The  use  of  analytical  coor¬ 
dinate  transformations  makes  the  method  efficient  and  the  results 
accurate.  Presently,  the  additional  pair  of  secondary  vortices 
indicative  of  Dean's  instability  for  curved  rectangular  ducts 
have  not  been  identified  for  the  polar  case.  A  systematic 
study  should  be  performed  for  the  fully  developed  flew  in  curved 
polar  ducts  using  the  appropriate  asymptotic  equations  to 
analyze  the  possibility  of  this  phenomenon  for  the  present 
configuration.  These  asymptotic  solutions  will  also  serve  as 
a  check  on  the  accuracy  of  the  fully  developed  flow  generated 
by  the  present  analysis.  In  fact,  these  asymptotic  results 
would  comprise  the  only  available  check  on  the  present  results 
since,  to  the  best  knowledge  of  the  present  investigators,  flow 
through  curved  polar  ducts  has  not  been  previously  analyzed. 

The  present  study  provides  the  first  set  of  available  results 
for  this  important  configuration.  To  render  this  research  even 
more  physically  meaningful,  results  are  to  be  generated  for  cases 
where  the  outer  wall  r  *  R-,  rotates  at  constant  angular  velocitv. 
Finally,  a  rotation  of  the  entire  duct,  with  this  outer  wall 
rotating  in  the  opposite  direction,  should  also  be  considered. 
This  configuration  very  closely  approaches  the  rotating  blade 
passages  in  turbcmachinery  applications. 


A  technical  abstract  (Ref.  10)  summarizing  the  results  of 
this  study  has  been  accepted  for  presentation  at  the  AIAA  Fluid 
and  Plasma  Dynamics  Conference  to  be  held  in  July  1979  in 
Williamsburg,  Virginia;  a  written  version  of  the  paper  is  to 
appear  later  as  an  AIAA  Preprint. 

1 . 3  Analysis  of  Laminar  Streamwise  Separation  via  a  Semi-Elliptic 

Mathematical  Model 

Sufficient  evidence  has  now  been  established  to  indicate 
that  incompressible  flow  seoaration  can  be  represented  by  mathe¬ 
matical  models  that  are  simpler  than  the  full  three-dimensional 
elliptic  Navier-Stokes  equations.  In  the  limit  of  infinite 
Reynolds  number,  the  triple-deck  theory  of  Stewartson  (Ref.  11) 
provides  an  asymptotic  analysis  for  separation  bubbles  of  small 
normal  extent.  For  finite  Reynolds  numbers,  the  various  simplified 
mathematical  models  considered  in  Refs.  12-14  lead  to  solutions 
that  are  in  good  agreement  with  the  corresponding  solutions  of 
the  full  Navier-Stokes  equations.  However,  all  these  finite- 
Re  analyses  have  been  in  the  framework  of  the  formulation 

for  two-dimensional  flows.  In  view  of  the  growing  interest  in 
three-dimensional  flow  and  the  difficulties  associated  with  the 
vorticity-stream  function  formulation  for  three-dimensional  flows, 
the  present  study  is  aimed  at  developing  r>  simplified  model  and 
an  efficient  solution  procedure  for  separated  flow  using  a 
velocity-pressure  formulation.  Although  results  are  oresently 
obtained  for  a  two-dimensional  model  problem,  the  entire  philosophy 
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in  the  study  is  directed  toward  a  three-dimensional  flow  analysis. 
Moreover,  for  this  two-dimensional  case,  results  were  also 
obtained  using  the  complete  Navier-Stokes  equations  in  order  to 
provide  an  accurate  evaluation  of  the  solutions  resulting  from 
the  simplified  model. 

The  model  problem  used  in  the  present  study  is  the  flow  in 
a  doubly  infinite  channel  with  an  asvmmetric  distortion  as 
shown  in  Fig.  12.  Surface-oriented  coordinates  for  this  channel 
are  the  conformal  coordinates  obtained  via  a  Joukowski  trans¬ 
formation.  The  straight  channel  with  parallel  walls  is  obtained 
as  a  special  case  of  this  configuration  by  the  variation  of  a 
parameter  in  the  transformation.  Clustering  of  the  coordinate 
curves  near  the  channel  walls  as  well  as  in  the  regions  of 
maximum  curvature  near  the  lower  wall  is  also  achieved  by 
analytical  coordinate-stretching  transformations.  For  certain 
configurations,  three  regions  of  separated  flow  are  anticipated 
to  occur  for  such  a  geometry;  these  are  indicated  as  regions  1, 

2  and  3  in  the  figure. 

The  flow  is  represented  mathematically  in  terms  of  the 
primitive  variables  (v,w,p)  and  the  conformal  coordinates  shown 
in  Fig.  12.  The  momentum  equations  governing  the  velocity  com¬ 
ponents  v  and  w  along  the  n  and  ;  directions,  respectively, 
are  parabolized  with  respect  to  the  streamwise  direction  ;. 

The  governing  equation  for  pressure  is  derived  by  forming  the 
divergence  of  the  momentum  equations,  and  using  the  continuity 
equation  indirectly  in  it.  Streamwise  ellipticity  is  maintained 
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in  the  pressure  equation;  hence,  the  tern  semi-elliptic  for  the 
proposed  model.  The  boundary  conditions  for  the  momentum  equa¬ 
tions  are  quite  simple  and  consist  of  zero  slip  at  the  walls 
and  zero  streamwise  variation  of  velocities  at  upstream  and 
downstream  infinity.  However,  any  realistic  velocity  profile 
may  be  prescribed  at  the  upstream  boundary.  The  boundary  condi¬ 
tions  for  the  pressure  equation  are  not  so  straightforward. 

They  must  be  derived  from  the  momentum  equations  evaluated  at 
the  boundaries  while  strictly  enforcing  zero  dilation  at  all 
boundaries.  Evaluation  of  the  normal  derivatives  appearing  in 
the  equation  at  the  boundaries  requires  special  care  as  it  has 
a  non-negligible  influence  on  the  solution  in  the  interior. 

The  initial  conditions  for  the  problem  consist  of  the  solution 
of  the  parabolized  moment-urn  equations  with  a  uniform  streamwise- 
pressure  gradient  corresponding  to  fully  developed  straight- 
channel  flow. 

The  finite-difference  equations  resulting  from  the  discre¬ 
tization  of  the  quasi linear ized  n-  and  ;-momentum  equations  on 
curves  of  ;  =  constant  are  solved  simultaneously  using  a 
generalization  of  the  Thomas  algorithm;  this  has  been  found  to 
be  more  efficient  than  conventional  block- tridiagor.al  solvers. 
Starting  with  prescribed  inlet  velocity  profiles,  the  streamwise 
marching  of  the  solution  proceeds  through  the  entire  computational 
field  up  to  the  downstream  boundary.  Thereafter,  the  pressure 
field  is  updated  by  a  prescribed  number  of  iterations  using  the 
Successive-Over  Relaxation  method  or  the  Alternating-Direction 
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Implicit  method.  The  local  dilation  as  well  as  the  mass  flow 
across  each  cross  section  of  the  channel  are  monitored  as  a 
check  on  the  accuracy  of  the  solutions.  Some  additional  remarks 
about  the  pressure  solution  are  worth  mentioning  here.  For  flow 
through  an  infinitely  long  channel,  the  pressure  becomes  unbounded 
as  the  downstream  boundary  is  approached.  Moreover,  since  the 
Poisson  equation  for  pressure  is  subject  to  Neumann  boundary 
conditions,  its  successful  solution  requires  the  satisfaction  of 
an  integral  relation  between  the  source  term  and  the  boundary 
values.  Numerical  consideration  of  this  relation  becomes 
impractical  for  an  unbounded  domain  of  integration.  To  circumvent 
these  difficulties,  the  problem  for  pressure  is  formulated  in 
terms  of  a  new  variable  that  represents  the  difference  between 
the  actual  pressure  distribution  and  the  fully-developed  straight- 
channel  pressure  distribution.  This  new  variable  remains  bounded 
everj'vhere;  moreover,  it  can  be  set  to  zero  at  the  upstream 
boundary,  thus  doing  away  with  the  Neumann  boundary  value  problem 
and  the  integral  relation  over  an  unbounded  domain. 

Results  were  first  obtained  for  the  degenerate  case  of  flow 
through  a  straight  channel  with  fully  developed  flow  conditions 
at  entrance.  The  final  converged  solution  was  obtained  in  a  single 
iteration.  Moreover,  forcing  several  additional  iterations 
showed  imperceptible  changes  in  the  solution,  thus  indicating  a 
stable  and  convergent  solution  procedure.  Next,  the  flow  in  a 
straight  channel  with  a  uniform  entrance  velocity  profile  was 
computed.  This  case  was  also  computed  using  the  fullv  elliptic 
Navier-Stokes  equations.  The  two  solutions  were  found  to  agree 
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very  well,  thus  conf inning  the  validity  and  accuracy  of  the 
semi-elliptic  model  proposed  in  the  present  study.  Calculations 
were  then  made  for  channels  with  various  degrees  of  distortion. 

As  the  distortion  increased,  reversed  flow  was  observed  in 
region  1  (Fig.  12)  of  the  channel.  Fine-mesh  calculations  also 
indicate  the  existence  of  a  very  small  senaration  bubble  in 
region  2.  In  all  cases,  the  results  of  the  semi-elliptic  model 
are  in  excellent  agreement  with  those  of  the  full  Navier-Stokes 
solutions.  Figure  13  shows  the  streamline  pattern  obtained  from 
the  computed  velocity-pressure  solution  for  a  typical  channel 
encountering  separation. 

Thus,  the  study  is  believed  to  have  served  its  primary 
objective  of  formulating  a  semi-e lliptic  model  and  an  accurate 
solution  procedure  for  laminar  incompressible  separated  flow. 
However,  these  results  are  considered  only  preliminary;  several 
points  are  being  presently  studied  further.  So  far,  no  particular 
effort  has  been  directed  toward  the  efficiency  cf  the  solutions. 

It  appears  that  a  more  consistent  finite-difference  formulation 
of  the  velocity-pressure  system  is  possible  with  use  of  a  stacgered 
grid  arrangement  than  with  the  regular  grid.  Therefore,  a 
staggered-grid  formulation  is  expected  to  increase  the  accuracv 
as  well  as  the  efficiency  of  the  solution.  Significant  gain  in 
solution  accuracy  is  also  expected  by  using  the  simplified  spline 
procedure  described  in  Ref.  15.  Hicher-order  accurate  fine-mesh 
solutions  will  enable  comparison  of  the  present  work  with  the 
asymptotic  solutions  obtained  by  Smith  (Ref.  16)  for  channel 
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flows.  The  true  test  of  this  research  lies  in  the  analysis  of 
three-dimensional  duct  flow  with  streanwise  separation.  An 
appropriate  model  problem  for  this  purpose  can  be  obtained  by 
placing  bounding  walls  a  finite  distance  aoart  across  the  width 
of  the  channel  in  Fig.  12. 

An  abstract  (Ref.  1?)  summarizing  this  work  has  been 
accepted  for  presentation  at  the  Symposium  on  Computers  in  Aero¬ 
dynamics  to  be  held  at  the  Polytechnic  Institute  of  New  York 
in  June  1979.  A  full  length  DaDer  is  nresentlv  under  preparation 
and  is  planned  to  be  submitted  for  publication  in  the  Inter¬ 
national  Journal  on  Computers  and  Fluids. 

1 . 4  Higher  Accuracy  of  Numerical  Solutions 

The  two-dimensional  Navier-Stokes  equations  for  incompressible 
flow  have  been  investigated  by  a  number  of  researchers.  The 
vorticity-strean  function  {w,v)  formulation  has  been  studied 
very  extensively  using  exdicit  as  well  as  implicit  numerical 
schemes  of  second-order  or  fourth-order  accuracy.  On  the  other 
hand,  the  velocity-pressure  (u,v,p)  formulation  had  been  solved 
mainly  by  explicit  numerical  procedures  so  far.  Recently,  implicit 
solution  procedures  have  also  been  apclied  to  the  (u,v,o)  formu¬ 
lation  (e.g.,  Ref.  6)  with  efficiency  and  accuracv  comparable 
to  or  better  than  those  for  the  (w,M  ^emulation .  The  orimarv 
difficulty  with  the  (u,v,p)  formulation  lies  in  the  determination 
of  a  pressure  field  that  can  lead  to  a  divergence-free  velocity 
field.  But,  it  is  the  (u,v,o)  formulation,  rather  than  the  (*,,.) 
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formulation,  that  admits  more  readily  to  extension  to  flows 
with  three-dimensionality,  turbulence,  etc.  With  a  trend  towards 
analysis  of  flow  problems  of  increasing  complexities,  the  accuracy 
and  efficiency  of  numerical  solutions  have  become  exceedingly 
important.  It  has  been  shown  (Refs.  18,  19)  that  the  use  of 
spline  approximations  in  the  discretization  of  the  (u,j)  dif¬ 
ferential  equations  has  several  important  advantages  over  con¬ 
ventional  finite-difference  approximations.  Derivative  boundary 
conditions,  irregular  boundaries  and  nonunifcrm  mesh  are  handled 
better  and  extension  to  higher-order  accurate  formulations  is 
more  straightforward.  Even  for  the  second-order  accurate  form, 
the  truncation  error  of  a  spline  formulation  is  lower  than  that 
for  other  second-order  schemes. 

The  principal  goal  of  the  present  study  is  to  employ  spline 
polynomials  in  the  numerical  solution  of  the  (u,v,p)  formulation 
of  the  Navier-Stokes  equations  and  to  analyze  the  accuracy  and 
the  efficiency  of  the  solutions.  The  incompressible  laminar 
viscous  driven  flow  in  a  square  cavity  is  selected  as  the  model 
problem  for  this  purpose.  This  is  a  true  Vavier-Stokes  problem 
with  no  admissible  boundary-layer-like  or  other  simplifying 
characteristics.  Furthermore,  experimental  as  well  as  accurate 
numerical  results  are  available  for  this  flow  problem. 

The  governing  equations  for  this  problem  are  the  momentum 
and  mass  conservation  ecruations;  only  the  convective  form  of  the 
equations  is  presently  used.  Since  the  pressure  does  not  appear 
as  a  dominant  variable  in  any  of  these  equations,  a  Poisson 
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equation  is  formed  for  the  pressure  using  the  momentum  and  con¬ 
tinuity  equations.  The  boundary  conditions  for  the  pressure 
equation  consist  of  normal  derivatives  of  pressure  determined 
from  the  momentum  equation  evaluated  at  the  boundaries.  As  in 
Ref.  6,  the  source  term  and  the  boundary  values  in  this  Neumann 
boundary-value  problem  are  required  to  satisfy  an  integral 
relation  necessary  for  a  convergent  solution.  Also,  the  pressure 
equation  is  formulated  directly  from  the  discretized  form  of  the 
momentum  equations  using  the  various  spline  relations. 

For  the  purpose  of  efficient  numerical  solution,  the  non¬ 
linear  terms  in  the  momentum  equations  are  quasi linearized  and 
the  resulting  equations  together  with  the  necessary  spline 
relations  are  solved  simultaneously  using  an  alternating-direction 
implicit  scheme  and  a  generalization  of  the  Thomas  algorithm 
(Ref.  20).  This  approach  is  computationally  much  faster  than 
conventional  (4  x  4)  block-tr idiagonal  solvers  and  is  stable  over 
a  wider  range  of  the  time  step  it.  Even  when  optimal  time  steps 
are  employed  in  the  momentum  as  well  as  the  pressure  equations, 
the  time-indeDendent  nature  of  the  continuity  equation  suggests 
that  more  than  one  pressure  iteration  be  performed  for  each 
moment-urn  iteration.  The  number  of  pressure  iterations  used 
varies  with  the  grid  size  emoloyed. 

The  spline  formulation  of  the  pressure  equation  oresents 
3ome  unique  difficulties  not  experienced  with  the  vorticity 
equation.  Even  though  both  i  and  p  and  their  derivatives  are 
singular  at  the  upper  corners  of  the  cavity,  the  ..-equation  is 
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not  required  at  these  singular  points,  whereas  the  pressure 
equation  is.  Various  formulations  have  been  attempted  for  the 
pressure  derivatives  at  these  singular  points,  including  a  cubic 
extrapolation  from  the  neighboring  points.  Considerable  improve¬ 
ment  is  possible  in  this  area  by  asymptotic  analysis  of  the 
corner  singularities  as  explained  later  in  this  section.  Never¬ 
theless  the  approach  used  presently  yields  reasonably  accurate 
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results,  with  the  velocity  divergence  in  the  order  of  10 

Results  have  been  obtained  for  several  values  of  Re  uo  to 
1000  and  compared  with  other  available  finite-difference  solutions 
as  well  as  spline  solutions  of  the  (^,v)  formulation.  As  indicated 
in  Ref.  19  for  (^,v)  spline  solutions,  the  present  (u,v,p)  spline 
solutions,  with  the  grid  sized  used,  are  more  accurate  than 
corresponding  finite-difference  solutions  even  with  a  finer  grid. 
For  purpose  of  an  accurate  comparison,  the  (w,.)  spline  solutions 
were  actually  computed,  by  the  present  authors  rather  than 
obtained  directly  from  Ref.  18.  Having  the  complete  (w,i) 
solution  available  also  enables  a  more  detailed  and  comprehensive 
comparison . 

Since  the  behavior  of  the  flew  variables  for  this  oroblem 
is  well  known ,  the  results  presented  here  are  only  those  corre¬ 
sponding  to  the  analysis  of  the  accuracy  of  the  solutions. 

Figure  14  shows  the  variation  of  vorticitv  _  at  the  midcoint 
of  the  moving  wall  with  variation  in  mesh  size  lx  (**  iy)  .  All 
three  sets  of  results  presented  indicate  second-order  accuracy. 

The  (u,v,p)  spline  solution  agrees  well  with  the  (w,^)  spline 
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solution,  and  both  these  solutions  have  much  smaller  error 
than  the  finite-difference  solution  obtained  by  the  factored 
ADI  scheme  of  Ref.  21. 

Figures  15  and  16  show  the  effect  of  mesh  size  on  the  stream 

function  .  at  the  vertex  center  of  the  cavitv  and  the  minimum 
max 

velocitv  u  at  the  centerline  of  the  cavitv.  The  (w,]/)  and 
*  mm 

the  (u,v,p)  soline  formulations  lead  to  j  and  u  that 

max  min 

agree  only  at  the  limit  of  vanishing  step  si2e.  In  one  formu¬ 
lation,  the  limiting  value  is  approached  from  above  while,  in 
the  other,  it  is  approached  from  below.  The  pressure  computed 
in  the  (u,v,p)  spline  formulation  is  also  indicated  to  be 
second-order  accurate  (Fig.  17).  Also  included  here  are  typical 
contour  plots  (Fig.  18)  for  the  cavitv  flow  at  Re  »  1000,  showing 
that  the  (u,v,p)  spline  formulation  developed  in  the  present 
research  is  capable  of  computing  this  high  Re  case  correctly 
even  with  a  uniform  (33  x  33)  grid. 

It  is  important  to  point  out  that  although  the  behavior  of 

2 

the  solution  appears  to  be  linear  with  respect  to  L  in  the 

figures  shown,  an  actual  calculation  of  this  exponent  shows 

that  the  exponent  is  frequently  far  from  2.  This  observation 

has  also  been  confirmed  for  several  ether  available  solutions 

that  graphically  indicated  second-order  accuracy.  The  authors 

believe  that  the  deoarture  from  exact  second-order  accuracy  is 

primarily  due  to  the  singularities  that  exist  at  the  two  upper 

corners  of  the  cavity.  Therefore,  an  effort  has  been  made  to 

2 

consider  nonuniform  motion  of  the  upper  wall  (e.g.,  U  3  sm  _x) 


in  order  to  eliminate  the  singularities.  More  importantly,  an 
effort  has  also  been  made  to  reformulate  the  dependent  variables 
of  the  original  cavitv-flow  problem  such  that  the  new  dependent 
variables  are  free  from  the  singularity.  The  local  singular 
solution  is  the  Stokes-flow  solution  determined  analytically 
as  an  asymptotic  series  (Ref.  22).  This  solution  vanishes  far 
from  the  corner  points  and  is  therefore  subtracted  from  the 
original  solution  to  yield  the  new  dependent  variables  which 
are  regular  everywhere.  Results  are  presently  being  comDuted 
using  these  new  variables.  These  results  are  expected  to  be  the 
most  accurate  and  efficient  results  possible. 

A  fourth-order  accurate  spline  formulation  of  the  (u,v,p) 
system  has  also  been  attempted.  Presently,  the  source  term  in  the 
pressure  equation  is  not  exactly  fourth-order  accurate.  The 
difficulty  arises  primarily  because  the  formulation  requires 
pressure  derivatives  at  the  singular  corner  points.  The  intro¬ 
duction  of  non-singular  dependent  variables  is  expected  to 
relieve  this  difficulty  also. 

A  technical  abstract  (Ref.  23)  based  on  the  results  of  this 
study  has  been  accepted  for  presentation  at  the  AIAA  Computa¬ 
tional  Fluid  Dynamics  Conference  to  be  held  in  Williamsburg, 
Virginia,  in  July  1979;  a  written  version  of  the  work  will  appear 
later  as  an  AIAA  paper. 

1 . 5  Efficient  numerical  Algorithms 

The  velocity-pressure  (u,v,p)  and  the  vorticitv-stream 
function  (.,;)  formulations  of  the  two-dimensional  unsteady 
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incompressible  Navier-Stokes  equations  have  been  solved  by  com¬ 
bining  two  very  efficient  and  accurate  numerical  techniques. 

The  velocity  transport  equations  in  the  primitive-variable 
(u,v,p)  formulation  and  the  vorticity  transport  equation  in 
the  vorticity-stream  function  formulation  are  solved  by  employing 
a  modified  form  (Ref.  6)  of  the  alternating-direction  implicit 
(ADI)  method  of  Peaceman  and  Rachford.  In  addition  to  these 
transport  equations,  both  formulations  contain  an  elliptic 
Poisson  equation  with  either  Dirichlet  or  Neumann  boundary 
conditions.  By  far,  the  greatest  computational  effort  is 
expended  in  obtaining  numerical  solutions  of  these  Poisson 
equations.  This  is  because  these  equations  are  truly  elliptic 
and  their  solution  time  increases  rapidly  with  the  number  of 
computational  points  employed.  For  3-D  unsteady  viscous  flow  at 
high  Reynolds  number,  reliable  numerical  solutions  must  employ  a 
reasonably  fine  computational  mesh.  Hence,  it  is  necessary 
to  develop  very  efficient  algorithms  for  solution  of  these 
Poisson  equations.  The  Hockney  (Ref.  24'  and  especially  the 
3uneman  (Refs.  25,  26)  Direct  Poisson  solvers  represent  two  such 
highly  efficient  solution  techniques  for  the  standard  "five 
point  star"  discretization  of  Poisson  equations  with  either 
Dirichlet  or  Neumann  boundary  conditions. 

The  nonlinear  convection  derivatives  occurring  in  the 
velocity  transport  equations  are  treated  by  quasi-linearization 
in  time,  thus  retaining  formal  second-order  space-time  accuracy 
without  the  necessity  of  iterating  the  nonlinear  terms.  However, 
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the  nonlinear  terms  in  the  vorticity  transport  equation  are 
merely  linearized  through  a  time  delay  of  the  stream  function, 
thus  resulting  in  a  difference  equation  which  possesses  second- 
order  spatial  accuracy  but  only  first-order  temporal  accuracy. 

Here,  the  temporal  second-order  accuracy  is  achieved  by  iterations. 
Also,  both  conservative  and  convective  finite-difference  forms 
of  the  transport  equations  are  investigated  with  regards  to 
computational  efficiency  and  the  accuracy  of  the  numerical 
solutions  obtained. 

Thus,  the  present  study  consists  of  a  detailed  numerical 
investigation  of  the  solution  of  a  model  two-dimensional  laminar 
incompressible  viscous  flow  configuration  by  a  total  of  eight 
numerical  procedures  (see  Table  1)  using  finite-difference  mesh 
spacings  of  (17  x  17)  to  (65  x  65)  grid  points,  over  a  range  of 
Reynolds  numbers  from  100  to  1000.  For  a  meaningful  comparison 
of  these  various  procedures,  Roache  (Ref.  27)  has  pointed  cut 
the  essential  requirement  that  the  common  testing  ground  be  a 
truly  Navier-Stokes  model  problem,  that  is,  one  without  any 
admissible  boundary-layer  or  other  such  simplifying  characteristics. 
Hence,  it  was  decided  to  investigate  the  problem  of  shear-driven 
flow  within  a  twc-dimens  lonal  scuare  cavity  as  the  model  croblem 
for  this  study  also.  The  corner  singularities  present  in  this 
flow  configuration  do  present  some  theoretical  as  well  as 
numerical  difficulties,  in  that,  the  numerically  estimated 
accuracy  of  the  results  reported  herein  is  degraded  below  second 
order  as  one  approaches  either  of  the  singularities.  To 
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circumvent  this  difficulty,  some  effort  has  again  been  made  to 
obtain  solutions  using  a  modified  velocity  distribution  at 
the  moving-wall  boundary,  namely,  u  =  sin^'x,  as  this  case  is 
expected  to  be  free  of  singularities. 

Numerical  solutions  for  the  vorticity,  stream  function, 
velocity,  and  pressure  fields  have  been  obtained  for  the  cavity 
flow  starting  from  a  no-flow  condition  and  proceeding  to  the 
fully  developed  recirculating  steady-state  flow  as  illustrated 
in  Fig.  19.  Here,  t  represents  a  nondimensional  unit  of  time. 

To  facilitate  the  comparison  between  solutions  using  the  vorticity- 
stream  function  (w,\H  variables  and  those  utilizing  the  velocity- 
pressure  (u,v,p)  variables,  a  regular  grid  system  is  presently 
employed.  Typical  results  for  the  axial  velocity  profile 
through  the  vertical  centerline  of  the  cavity  have  been  plotted 
in  Figs.  20  and  21.  As  seen  in  Fig.  20,  the  results  of  the  present 
study  using  the  (~,v)  system,  with  uniform  (33  x  33)  and  (65  x  65) 
grids,  are  in  good  agreement  with  the  results  available  in  the 
literature.  In  Fig.  21,  the  present  (u , v , p) -system  results  are 
compared  with  the  (w, v) -system  solutions.  The  slight  discrepancy 
seen  in  this  figure  seems  to  be  a  direct  consequence  of  the  coarse 
grid  used  fcr  the  (u,v,p)  system.  It  is  important  to  note  that 
the  current  convective  ADI-Buneman  (w,t)  solution  represents 
an  order  of  magnitude  improvement  in  computational  efficiency 
over  previously  published  results  (see  Ref.  28). 

A  heuristic  accuracy  study  following  the  lead  of  Blottner 
(Ref.  29)  and  illustrated  in  Fig.  22,  indicates  that  the  formal 
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second-order  spatial  accuracy  of  these  Navier-Stokes  solvers 
degrades  somewhat  as  one  approaches  the  locations  of  the  corner 
singularities  of  this  model  problem.  This  heuristic  accuracy 
analysis  uses  Richardson's  extrapolation  method  and  the  results 
for  the  (17  x  17),  (33  x  33),  and  (65  x  65)  grids  to  estimate 

the  infinitesimal-grid  solution.  This  infinitesimal-grid 
estimate  is  then  employed  to  yield  an  estimate  of  the  percentage 
error  level  existing  in  a  particular  solution.  A  log-log 
plot  of  the  absolute  values  of  the  oercentaqe  error  versus  number 
of  grid  points  along  the  x-axis  is  then  executed.  If  these 
points  lie  along  a  straight  line  whose  slope  is  -2,  then  the 
assumed  second-order  spatial  accuracy  is  confirmed;  deviation 
from  this  stra lght- 1  me  behavior  represents  degradation  in 
the  accuracy  level. 

In  summary,  efficient  and  accurate  numerical  solutions  of 
the  Navier-Stokes  equations  for  flow  in  a  driven  cavity  have 
been  obtained.  The  Poisson  equations  with  Dirichlet  and  Neumann 
conditions  in  the  (w,ii)  and  (u,v,p)  systems,  respectively,  have 
been  solved  using  the  fast  direct  methods  of  Buneman  as  well 
as  Hockney.  A  careful  comparative  study  of  the  numerical 
solutions  obtained  with  various  algorithms  is  carried  out; 
extensive  results  are  presently  being  generated.  A  couple  of 
problem  areas  have  been  identified  as  potential  candidates  for 
further  work;  these  are  briefly  outlined  here:  (i)  The  corner 
singularities  can  be  removed  by  employing  analytical  asymptotic 
solutions  at  the  corner.  Introduction  of  new  dependent  variables 
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will  not  only  improve  the  accuracy  but  will  aid  also  in  improving 
the  efficiency  of  the  solution.  Some  effort  has  already  been 
initiated  towards  this  analysis.  (ii)  The  (u ,v ,p) -system 
solutions  are  sought  with  a  view  that  this  system,  is  easily 


extendable  to  three-dimensional  flows.  However,  the  present 


A 

solutions  show  that  the  dilation  D  -  u 

x 


v  does  not  approach 

y 


a  desired  low  level.  It  now  appears  that  use  of  a  "staggered" 


grid  as  suggested  by  researchers  at  Los  Alamos  Laboratory  (Ref.  30) 


will  lead  to  desired  level  of  dilation  in  the  field.  Based  on 


the  study  carried  out  so  far,  three  different  types  of  formula¬ 
tions  for  the  discretized  problems  are  identified.  These  are 
associated  with  achieving  desired  low  levels  of  dilation  and 


are  referred  to  as  weakly,  mildly  and  strongly  solenoidal 
formulations;  these  are  currently  being  studied  further. 

An  abstract  (Ref.  31)  based  on  the  results  of  this  work 


has  been  submitted  for  presentation  at  the  Second  International 
Conference  on  Mathematical  Modelling  to  be  held  in  St.  Louis, 
Missouri,  in  July  1979;  a  written  version  of  the  paper  will 
then  appear  later  in  the  Conference  Proceedings. 
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2 .  FUTURE  EFFOR" 


During  the  remainder  of  the  contract  Deriod,  the  various 
projects  described  in  the  preceding  section  will  be  continued 
and  brought  to  a  state  of  completion.  This  will  include  the 
preparation  of  technical  reports  or  papers  describing  the 
research  in  reasonable  detail,  in  addition  to  completing  the 
research  itself.  For  example,  in  the  study  of  duct  flows, 
rotation  of  the  outer  wall  for  curved  polar  ducts  will  be 
considered  in  order  to  simulate  a  rotating  annular  cascade 
and  study  the  effect  of  this  wall  rotation  on  the  flow  in  polar 
ducts.  In  the  area  of  improvement  of  the  accuracy  of  numerical 
solutions  and  the  efficiency  of  solution  algorithms,  an  effort 
will  be  made  to  develop  an  asymptotic  solution  for  the  flow 
near  the  upper  corners  of  the  driven  cavity.  This  will  lead 
tc  the  formulation  of  new  dependent  variables  that  are  non¬ 
singular  at  these  corners.  Computation  of  the  flow  solution 
in  terms  of  these  singularity-free  variables  is  expected  to  lead 

to  a  more  accurate  and  efficient  numerical  solution.  For  the 

2 

same  purpose,  the  cavity-flow  configuration  with  U  *  sin  ttx 
will  also  be  analyzed  more  fully  as  this  configuration  does  not 
contain  any  singularities.  Finally,  a  staggered-grid  finite- 
difference  discretization  as  well  as  a  fourth-order  accurate 
spline  formulation  of  the  (u,v,o)  system  will  be  studied  further 
in  order  to  obtain  accurate  flow  solutions  of  the  cavity-flow 
problem  with  minimum  dilation  for  the  velocity  field. 


The  work  discussed  so  far  is  for  laminar  flow. 


However , 


most  flows  in  practical  applications  are  turbulent  in  character. 
Therefore,  a  concurrent  effort  has  been  also  initiated  for  the 
study  of  turbulent  flow.  This  effort  is  presently  at  a  rather 
preliminary  stace;  hence,  it  was  not  included  in  the  discussion 
in  the  preceding  section.  The  two-equation  turbulence  model 
of  Launder  and  Spalding  (Ref.  32)  has  been  used  by  Sokhey, 

K.  Ghia  and  U.  Ghia  (Ref.  33)  to  study  the  three-dimensional 
turbulent  flow  inside  curved  ducts  of  reoular  cross-sections, 
with  the  proper  wall  functions  for  the  regions  near  the  duct 
walls.  However,  the  accuracy  as  well  as  the  efficiency  of  these 
numerical  solutions  needs  to  be  studied  in  detail.  Analytical 
transformations  of  the  form  used  in  Ref.  6  will  be  used  to 
provide  the  desired  clustering  of  grid  points  in  the  regions  of 
high  gradients  near  the  walls.  Solution  of  the  flow  equations 
and,  in  particular,  of  the  Poisson  equation  for  pressure  with 
Neumann  boundary  conditions,  in  these  transformed  variables 
is  not  a  trivial  task.  The  results  of  the  turbulent  flow  cases 
will  be  verified  with  the  experimental  data  in  order  to  demon¬ 
strate,  in  a  convincing  manner,  the  usefulness  of  the  turbulence 
model  and  the  associated  analysis  developed. 

The  wall  regions  will  also  be  treated  in  a  unified  manner 
by  employing  the  Jones  and  Launder  turbulence  model  (Ref.  34) 
for  these  low-Reynolds  number  regions.  The  results  obtained  by 
the  wall-function  method  and  the  low-Re  modelling  will  be  com¬ 
pared  in  order  to  evaluate  these  two  approaches  for  treating 
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the  wall  regions.  Finally,  an  effort  will  also  be  made  to 
introduce  compressibility  in  the  analysis.  For  this  purpose, 
the  two-equation  turbulence  model  of  Launder  and  Spalding 
(Ref.  32)  will  be  used  to  study  a  typical  compressible  flow 
problem,  such  as  the  flow  exhausting  from  a  rectangular  nozzle 
for  which  the  corresponding  laminar-f low  analysis  and  solution 
procedure  are  already  available. 
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PERSONNEL 


3  . 

The  research  reviewed  in  this  annual  report  represents 
the  efforts  of  the  two  principal  investigators  for  task  II 
of  the  subject  research  project:  Drs .  K.N.  Ghia  and  U.  Ghia. 
Also,  Mr.  G . A .  Osswald  and  Mr.  C.'T>.  Shin,  graduate  research 
assistants  in  the  Aerospace  Engineering  Department,  have 
assisted  in  the  research.  Mr.  Osswald's  efforts  are  directed 
primarily  in  the  area  of  developing  and  employing  direct 
Poisson  solvers  in  order  to  improve  the  efficiency  of  numerical 
solutions  of  the  Navier-Stokes  equations.  Mr.  Shin's  contri¬ 
butions  have  been  in  the  use  of  spline  approximations  in  the 
aescretization  of  the  velocity-pressure  formulation  of  the 
Navier-Stokes  equations.  These  efforts  of  Mr.  Osswald  and 
Mr.  Shin  are  leading  to  their  Master's  degree  theses  which  are 
expected  to  be  completed  during  the  remainder  of  the  contract 
period . 
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FIGURE  1 .  SURFACE-ORIENTED  COORDINATES  FOR  VISCOUS  FLOW  THROUGH  A 
TYPICAL  ROTOR  CASCADE. 
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Aspect  Ratio  ®  1.0,  Area  Ratio  ®  0.8.  tmax  =  0.0*28,  Re  *  028 


(z/Re  D) 

AT  ION  OF  CENTERLINE  VELOCITY,  MAXIMUM  VELOCITY  AND  MEAN-VISCOUS  PRESSURE  DROP 
NORMALIZED  AXIAL  DISTANCE  FOR  POLAR  JOUKOWSKI  DUCT. 
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FIGURE  8  .  COMPARISON  OF  STREAM,  WISE  VARIATION  OF  CENTER 
LINE  AXIAL  VELOCITY  AND  MEAN  VISCOUS  PRESSURE 
FOR  CURVED. POLAR  AND  CURVED  SQUARE  DUCT 
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FIG.  9.  EFFECT  Or  DEAN  NUMBER  ON  SKIN-FRICTION  IN 

STREAMW 1 SE  DIRECTION  IN  A  CURVED  POLAR  DUCT. 
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FIGURE  10b.  DEVELOPMENT  OF  AXIAL  VELOCITY  PROFILES  IN  A  CURVED  POLAR  DUCT 
K  =  630,  R  =  2520,  R  =  16,  AR  -  2/3 


FIGURE  11a.  FULLY  DEVELOPED  SECONDARY  VELOCITY  PROFILES 
FOR  A  CURVED  POLAR  DUCT, 

K  =  100,  R/D  =  S3. 3,  R  =  913,  AR  =  2/3 
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D  SECONDARY  FLOW  VELOCITY  PROFILES 
i_AR  DUCT  K  *  63C,  R  =  16,  AR  =  2/3 


FIG.  15.  STREAMLINES  FOR  SEPARATED  FLOW  IN  CONSTRICTED  CHANNEL. 
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i.  15.  ACCURACY  OF  SPUNE-2  SOLUTIONS  AND  FINITE-DIFFERENCE 
SOLUTIONS  FOR  STREAM  FUNCTION  AT  VORTEX  CENTER 
OF  DRIVEN  CAVITY. 
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/•)  SPLINE  AND  (u,v,p)  SPLINE  SOLUTIONS 
VITY  CENTERLINE. 


HORIZONTAL  VELOCITY  PROFILE  THROUGH  CAVITY  CENTER 


G.  22.  ACCURACY  STUDY  p0P  ( <*• /. ) -SVSTEM  SOLUTIONS 


FOR  DRIVEN-CAVITY  FLOW. 


